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Abstract 

Quantum systems with sublevel structures prevent full population inversion from one manifold of 
sublevels to the other using strong ultrafast resonant pulses. In this work we explain the mechanism 
by which this population transfer is blocked. We then develop a novel concept of geometric control, 
assuming full or partial coherent manipulation within the manifolds and show that by preparing 
specific coherent superpositions in the initial manifold, full population inversion or full population 
blockade, i.e laser-induced transparency, can be achieved. In particular, by parallel population 
transfer we show how population inversion between the manifolds can be obtained with minimal 
pulse area. As the number of sublevels increases, population inversion can overcome the pulse area 
theorem at the expense of full control over the initial manifold of sublevels. 
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In this work we are concerned with intrinsic properties of the dynamics of systems with 
manifolds of sublevels, described by two (or more) quantum numbers, that hereafter will 
be generically called quantum structures. From the point of view of controlling the sys¬ 
tem dynamics, quantum structures pose several interesting problems. Quantum control[l] 
typically implies the ability to manipulate interfering pathways, which increases with the 
number of levels that participate in the dynamics as long as the system is controllable [2]. A 
multi-level structure would therefore offer more control opportunities at the expense of the 
ability to manipulate within the substructures. Our general goal is to investigate whether 
quantum structures limit, or conversely help, in controlling the system. In this paper we 
will be concerned with coarse-grained goals, where the objective of the control will be the 
state of the manifold given by the hrst quantum number, not the detailed state of the sub- 
levels. In hnding the best possible controls we will assume that the substructure is partially 
controllable, that is, that given some constraints, any possible wave function within a subset 
of the sublevels can be prepared [3] . Building on this assumption we will develop a geomet¬ 
ric control approach that allows hnding the optimal initial wave functions that maximize 
the yield of the desired process. This procedure does not prescribe an optimal held, but 
implicitly assumes that a held can be found, and makes full use of the quantum structures. 

Let us consider a simple and very general process in systems with a congested spectrum: 
absorption from the ground, initial manifold to the excited, target manifold, by means of 
a strong ultrashort pulse, with a bandwidth much larger than the energy spacing of the 
sublevels within the hne structure. As a molecular example, we can conceive controlling 
an electronic transition using a broadband pulse, where the goal is to invert the population 
to the excited state, regardless of the vibrational populations [4]. Although one may think 
that the substructure, specially in the case of very diherent associated time-scales, does not 
ahect the overall transition, in the strong held case the opposite occurs. The unpopulated 
levels of the hne structure induce Stark shifts [5] and create ehective detunings from the 
resonance that limit the extent of the Rabi oscillations. That is, assuming that all the 
diherent sublevels are dipole allowed, regardless of the strength of the pulse, the maximum 
population that can be reached is typically much smaller than one. 

A simple theoretical model explains this observation. Let us hrst assume that the diherent 
sublevels within each manifold are degenerate, AE Au ~ 0, and that all transient dipoles 
are equal. Then the equations of motion for every sublevel of the excited manifold \e,k), 
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and for every sublevel in the initial manifold \g,j), are the same. For a resonant transition: 
cij = iQ{t) ^k{t)/2 and bn = iflit) where fl{t) is the Rabi frequency. However, 

the initial conditions in the manifold \g,j), where we assume a single state, \g, 1), is initially 
populated, break the symmetry so that three different probability amplitudes describe the 
dynamics. We write, 

Nf Ni-l 

I'J'W) = ait)\9, 1) + Kt) |e, k) + c(t) Y \9^f) (1) 

k j' 

where the prime indicates that the initial state is excluded from the summation in the |i, j) 
manifold, and h{t) and c{t) are mean probability amplitudes, which in fact behave exactly as 
every sublevel amplitude. We dehne now the collective excited \E) and Raman \R) states 
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which together with the initial state |z) = \g, 1) form an orthonormal basis such that |'h(t)) = 
a{t)\i) + B{t)\E)+C{t)\R) with B{t) = Y^b{t), C{t) = y/Ni — lc(t), and \a{t)\‘^ + \B{t)\‘^ + 
|C(f)p = 1. The effective Hamiltonian in this basis is 
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H := (|^)(E| + y/W-l|^)(R| + c.c.) 


(4) 


with simple analytic eigenvalues and eigenvectors {aka dressed states) [6]. When W = 2 
the wave function dynamics is particularly interesting, with a population in the manifold of 
excited states given by 
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where 9{t) = jY^2{t')dt', such that ^(oo) is the pulse area, A. A maximum of 50% 
population can reach the excited state, whereas there is Rabi flopping (at twice the period 
of oscillation) between the initial state |i) (|i, 1)) and \R) (|i,2)). 


Pn{t) = \{RmtW = ^^^* 
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FIG. 1: Population histories for a system with (a) Ni = 2 {Nf = 1), (b) Ni = 7 {Nf = 7). In both 
cases the pulse area is ^ = Svr and the energy difference between the levels is AE = 0.4r“^. 

i.e. there is a very efficient Raman Stokes transition. Increasing the number of sublevels in 
the initial manifold only blocks the population transfer more efficiently. For large iV,, we 
obtain a population in the excited manifold of 

PE{t) = (7) 

Unlike in Ni = 2, for large iVj the maximum Raman populations depend also on N~^. 
Larger pulse areas only increase the frequency of the oscillations, but the initial state is 
mostly decoupled. The population dynamics behave as in an off-resonant excitation, with 
an effective detuning created by the Antler-Townes splitting[7] between the \E) and |i?) 
states, coupled by an y/Ni — 1 stronger Rabi frequency than the initial state with the \E) 
state. 

Fig.l shows the population dynamics for different sublevel structures, using time-scaled 
units. Case 1 refers to A, = 2, AE = 0.4r“^ < Au = 41n2 (where Au is the bandwidth 
of the Gaussian pulse, with scaled width, t = 1), Nf = 1 and the peak Rabi frequency is 
Go = 2. In case 2 we use Ni = Nf = 7 with the same energy splitting and pulse parameters 
as before. The results show (i) efficient Raman transfer for the first case and (ii) population 
locking in the second case, qualitatively in agreement with Eqs.(6) and (7) even for non- 
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degenerate structures. The main effect provoked by the energy splittings is to allow more 
population flow to the most excited sublevels of the initial and hnal manifolds, because 
the energy difference partially off-sets the effective detuning. But the effect is too small to 
qualitatively change the dynamics. 

Is it possible to optimize the pulse parameters to increase the efficiency of the population 
transfer? Clearly, as long as the initial state is a single sublevel Eq.(7) limits the maximum 
population that can be transferred using transformed-limited pulses. The control requires 
manipulation of the initial wave function. We will assume that the initial manifold can be 
manipulated before Q{t) acts such that we have full controllability within a given subset 
of states. Typically this requires the use of laser pulses of very different frequencies. In 
molecular physics, several control schemes have been proposed that imply creating coher¬ 
ences in the initial electronic state by means of infrared pulses before the optical held is 
used[8]. However, instead of explicitly hnding these pulses, in this work we will develop a 
geometrical approach. 

We want to maximize the population on the hnal manifold at time T, given by the 
functional J-, 

Nf 

7 = ^(d'(ti)|U(tj,T;H)|e,m)(e,m|U(T,ti;H)|T(ti)) (8) 

m 

for hxed H(t), with respect to changes in the initial wave function |\k(tj)), where is the 
initial time. This amounts to hnding U(tj, 0; ^i)\g, 1) = |\k(tj)). Instead of explicitly hnding 
the new held, we assume controllability and use a variational approach to simply obtain 
the rotation matrix R^l^f, 1) = \'ipi) = aij\g, j) (where the sum can be constrained to 
a subset of the levels of the initial manifold, i.e. < W) such that 7 is maximal. We 
therefore substitute |T(tj)) by |V’i) in Eq.(8). The optimization is purely geometrical and 
one can use the Rayleigh-Ritz approach. Thus we construct the matrix F with elements 

^7 = X^(bj|U(0,T;H)|/,m)(/,m|U(T,0;H)|i,fc) (9) 

m 

Restricting I'lpi) to be normalized (equivalently, Rj to be unitary) we obtain the secular 
equation, Fl-i/^j) = The solutions are the eigenvectors of F which give the yields of 

population transfer Xi- 

In Fig.2 we compare the results of the optimization with the yields obtained from the 
initial state |z) for degenerate structures AE = 0 and with constant energy spacing between 
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FIG. 2: Maximum population transfer to the excited manifold at final time as a function of pulse 
area for a system with (a) Ni = 2 and (b) Ni = 5. In both cases Nf = Ni. Dotted lines represent 
results for degenerate structures while solid lines give results for structures with energy spacings 
AE = 0.4r-b 

adjacent levels AE = 0.4r“^, for two different manifolds, = iVj = Nf = 2 and N^. = 
Ni = Nf = 5. The optimal yields for degenerate structures follow the pattern of Rabi 
oscillations coinciding with the population transfer from the initial state, but with full 
population transfer at multiples of tt of the pulse area. For non-degenerate structures the 
minima of these oscillations does not drop to zero but increases with the number of levels, 
despite the N~^ factor in Eq.(7). 

It is interesting to analyze generic features of the optimized initial states. When the 
extended pulse area Ae = ^^NiNfA is an odd multiple of tt, the optimal initial states have a 
very clear structure: All their coefficients are equal. This result can be explained analytically 
for the degenerate structure, following exactly the same steps as in Eq.(l). We now define 

Np _ Nf Nu 

|T(f)) = a{t) ^ \g, 1) + h{t) ^ |e, k) + c{t) ^ \gA) (10) 

j k j 

where a{t) is the mean amplitude of all Np initially populated levels in the initial manifold, 
whereas Nu = N^ — Np is the set of unoccupied states. Together with the previously defined 
|i?) and \E) collective states, the collective initial state 

. Np 

= ( 11 ) 
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forms a orthonormal set such that the Hamiltonian can be written as 


H := -1 (^V^|/)(E| + ^u\E){R\ + c.c.) (12) 


with the same eigenvalues as before. Given the wave function |^(t)) = A{t)\I) + B{t)\E) + 
C{t)\R) with A(t) = ^/Npa{t), the probability of reaching the excited manifold is 


PEit) 


NpN, 

+ Nl 


sin^ 



(13) 


Whenever ^e(^) is an odd multiple of tt, full population inversion can be achieved if all the 
sublevels of the initial manifold are equally populated and in phase. For different energy 
spacings other choices of phases and populations give better results, but the populations 
are always almost equal. On the other hand, it is simple to proof that when the initial 
probability amplitudes are all out of phase (such that Oi(0) = 0) then R is minimized 
and perfect transparency can be achieved, i.e. the population in the excited manifold is 
zero at all times[9]. Similar results are obtained with nondegenerate structures although 
the transparency is no longer perfect. Since there are many more possible solutions that 
minimize the yield (exactly Ni — 1 orthogonal eigenvectors for the degenerate structure) than 
those that maximize the yield (a single solution for the degenerate case) the set of eigenvalues 
hlls from below and the subspace of population transfer is of very small dimension. From 
the point of view of quantum controllability, population transfer is a difficult problem. 
Nevertheless, other yields greather than zero can be achieved when AE ^ 0. 

We now return to the original question: Do quantum structures help to control the dy¬ 
namics? For population inversion between two manifolds, the problem is easily controllable 
by means of proper Rabi oscillations in the simplest case, with W = Nf = 1. The existence 
of a substructure creates an effective detuning that reduces the oscillations. Therefore, the 
larger the energy spacing AE is, the smaller the back-effect of the unpopulated states and 
the more the system resembles a simple 2-level system. However, by manipulating the state 
within the initial substructure one can regain full Rabi oscillations and reach a regime where 
the pulse area theorem is overcome with almost full controllability regardless of the pulse 
area for very large W- Fig.3(a) shows how the optimized yield with respect to changes in the 
initial state practically achieves full population inversion whenever W = Ni is large. Fig.3(b) 
shows how the hrst (and second) minima of the optimal yield increase with W- Only in the 
degenerate case W does not play any role and the dynamics is less controllable. Moreover, 
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FIG. 3: Maximum population transfer to the excited manifold at final time as a function of the 
extended pulse area for systems with different number of sublevels Ni = Nf = N and fixed 
AE = 0.4t“^. In (a.) N = 2 (dotted line), N = 5 (green line), N = 20 (orange line). In (b) we 
show how the first and second minima of the optimized yield increases with N^. As shown in (a), 
the first minima for Ni = 20 occurs practically at the first maxima for smaller N. 

since the extended area increases with the number of sublevels, population inversion can be 
achieved with relatively weak pulses[10]. In the given example, with = Ni = Nf = 20, 
full inversion is obtained already with a pulse area of A ~ tt/IO (see Fig. 4). This is a 
consequence of parallel transfer. 

On the other hand, one should remember that the ability to optimize the yield with 
increasing Ni (for hxed AE) is at the expense of a hner optimization of the initial state, 
where the optimal solution is typically a single one, while the “robust” subspace of near¬ 
zero eigenvalues occupies practically all the space of solutions. What happens when the 
dimensionality of the subspace that is controlled is smaller than the space that is initially 
accesible or, in other words, how are the solutions deteriorated when W < W? In Fig.4 we 
show the yield of population transfer as a function of the pulse area for different “control 
subspaces”. Now W = 20 but the controller has only access to the hrst W sublevels (or 
to the odd numbered sublevels, due to e.g. a unspecihed selection rule or symmetry). The 
case Nc = 20 gives maximum yield while W = 1 implies no control over the initial state. 
The yields are deteriorated as the ability to control the system decreases and this effect 
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FIG. 4; Maximum population transfer to the excited manifold at final time as a function of the 
pulse area for a system with Ni = Nj = 20 sublevels with control on different subsets of the initial 
state, Nc- 

cannot be overcome by increasing the pulse area. On the contrary, often best results are 
often obtained with Ae = Moreover, adding external constraints, such as access to only 
odd number of levels, returns lower values of the yields. 

In summary, we have shown that quantum substructures may hamper the success of 
population transfer in multi-level systems. By engineering the initial state one can avoid 
the detrimental effects and partially correct the Rabi oscillations of the yield forced by the 
pulse area theorem. This is not achieved by brute force (increasing the pulse areas) but by 
preparing quantum superposition states that cancel the detrimental Raman transitions and 
lead to parallel transfer. Full population blockade and in fact laser induced transparency 
can also be achieved in similar manners. The control over the dynamics increases with 
the ability to manipulate every sublevel of the quantum substructure and is substantially 
reduced when there is limited control over the sublevels. 
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